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We calculate moments of the 0(af) heavy flavor contributions to the Wilson coefficients of the 
structure function /^(jc, Q 2 ) in the region Q 2 ^> m 2 . The massive Wilson coefficients are obtained 
as convolutions of massive operator matrix elements (OMEs) and the known light flavor Wilson 
coefficients. The calculation of moments of the massive OMEs involves a first independent re- 
calculation of moments of the fermionic contributions to all 3-loop anomalous dimensions of the 
unpolarized twist-2 local composite operators stemming from the light-cone expansion. The ex- 
pressions for the massive OMEs are now known for general values of the Mellin variable apart 
of one coefficient in the constant term, for which only a series of moments is computed. 
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1. Introduction 

Deep-inelastic scattering (DIS) of charged or neutral leptons off proton and deuteron targets, in 
the region of large enough values of the gauge boson virtuality Q 2 = —q, allows to measure the 
leading twist parton densities of the nucleon, the QCD-scale Aqcd, and the strong coupling con- 
stant a s {Q ) = a s (Q 2 )/(4n), to high precision. For unpolarized DIS via single photon exchange, 
the double-differential cross section can be expressed in terms of two inclusive structure functions 
Fi,l(x, Q 2 )- These decompose for twist z = 2 into a Mellin convolution of non-perturbative mass- 
less parton densities fj(x,n 2 ) and the perturbative Wilson coefficients ^ j^^ix.Q 2 j \x 2 ,rn 2 k l \x 2 ). 
The latter describe the hard scattering of the photon with a massless parton.They are given by the 
sum of the purely light - denoted by C y ( 2 ,l) ~~ an d heavy flavor contributions, H j nt) ■ Here k = c, b 
and j = q, g, depending on the type of process one considers, x denotes the Bjorken scaling vari- 
able. Especially in the region of smaller values of Bjorken-x, the structure functions contain large 
cc-contributions of up to 20-40 %, denoted by /^(x, <2 2 ). The perturbative heavy flavor Wilson 
coefficients corresponding to these structure functions are known at NLO semi-analytically in x- 
space [[p. Due to the size of the heavy flavor corrections, it is necessary to extend the description 
of these contributions to 0{a? s ), and thus to the same level which has been reached for the massless 
Wilson coefficients [||]. 

A calculation of these quantities in the whole kinematic range at NNLO seems to be out of 
reach at present. However, in the limit of large virtualities Q 2 , Q 2 > 10m 2 . in the case of F 2 cc (x, Q 2 ), 
one observes that F^Ax, Q 2 ) are very well described by their asymptotic expressions [Q] neglecting 
power corrections in m 2 / Q 2 , cf. also [Q]. In this kinematic range, one can calculate the heavy flavor 
Wilson coefficients analytically. This has been done for F 2 cc (x, Q 2 ) to 2-loop order in [||, |5|] and for 
F[ c (x, Q 2 ) to 3-loop order in [^J. Note that in the latter case, the asymptotic result becomes valid 
only at much higher values of Q 2 . The asymptotic expressions are obtained by a factorization of the 
heavy quark Wilson coefficients into a Mellin convolution of massive OMEs Aj^ and the massless 
Wilson coefficients Cjj, if one heavy quark flavor of mass m and n f light flavors are considered. In 
the present paper, we report on the calculation of the massive OMEs Ajk to 3-loop order for fixed 
even moments of the Mellin variable N, cf. [f7|] for details. We further calculate the OMEs which 
are required to define heavy quark parton densities in the variable flavor number scheme We 
also obtain moments of the terms oc T F of the 3-loop unpolarized anomalous dimensions yj. Our 
results agree with those obtained in [§]. Since the present calculation is completely independent by 
method, formalism, and codes, it provides a strong check on the previous results. 

2. Heavy Flavor Operator Matrix Elements 

The heavy flavor Wilson coefficients for a single massive quark may be expressed as 

H S.PS,Nsf Q 1 ™ 2 \_„S,?S ( Q 2 rr^_\ s,?S,Hs( nu 

M i,(2,L) ^ ^2 > ^2 ) ~ M j,(2,L) \ *> ^2^2) + L U2,L) \*> ^ > ^ ) » W 

where the photon couples to a light (L) or heavy (H) quark line, respectively. Further S stands 
for the flavor-singlet contributions, which are separated into a pure-singlet (PS) and non-singlet 
(NS) part via S = PS + NS. The factorization formula for the inclusive Wilson coefficients reads in 
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Mellin space, 



r^S,PS,NS,as f f. T „ Q 2 m2 \ V 4 S,PS,NSA. m2 \ n S,PS,NS f AT ^ Q 2 \ n ~ 



Here /I refers to the factorization scale between the heavy and light contributions in ^},(2,i) an d 
'as' denotes the limit Q 2 m 2 . The C ; - ( 2 ,l) ^ precisely the light Wilson coefficients and describe 
all the process dependence. The arguments («/), («/ + 1), indicate at how many light flavors the 
respective quantities have to be taken. This factorization is only valid if the heavy quark coefficient 
functions are defined in such a way that all radiative corrections containing heavy quark loops are 
included. Otherwise (2.2) would not show the correct asymptotic <2 2 -behavior The mass de- 
pendence is given by the process independent massive OMEs A, 7 , which are the flavor-decomposed 
twist-2 operator matrix elements 

2 oo 2 

» v^2' N ) = ^ °k l^ = 5 «+L a A ; {^2' N ) ■ (2 - 3) 

Here, i denotes the external on-shell particle (i = q,g) and Ok stands for the quarkonic (k = q) or 
gluonic (k = g) operator emerging in the light-cone expansion. The subscript H indicates that we 
require the presence of heavy quarks of one type with mass m. The logarithmic terms in m 2 / ji 2 are 
completely determined by renormalization and contain contributions of the anomalous dimensions 
of the twist-2 operators. Thus at NNLO the fermionic parts of the 3-loop anomalous dimensions 
calculated in Refs. [Q| appear. All pole terms of the unrenormalized results provide a check on our 
calculation and the single pole terms allow for a first independent calculation of the terms oc T F of 
the 3-loop anomalous dimensions. 

In case of the gluon operator, the contributing terms are denoted by A gq _Q and A gg Q. For the 
quark operator, one distinguishes whether the operator couples to a heavy or light quark. In the 
NS-case, the operator, by definition, couples to the light quark. Thus there is only one term, A^q. 
In the S and PS-case, two OMEs can be distinguished, {A^q, A^ q} and {A^, Ag g }, where, in 
the former case, the operator couples to a light quark and in the latter case to a heavy quark. 

Eq. (2.2) allows to calculate the heavy flavor Wilson coefficients in the limit Q 2 ^> m 2 up to 
0{o? s ) by combining the results obtained in Ref. [Q] for the light flavor Wilson coefficients with the 
3-loop massive OMEs which are computed in this work 

A related application of the heavy OMEs is given when using a variable flavor number scheme 
to describe parton densities including massive quarks. The OMEs are then the transition functions 
going from tif to n/ + 1 flavors. One thus may define parton densities for massive quarks, see 
e.g. Ref. [§]. This is of particular interest for heavy quark induced processes at the LHC, such as 
cs — > W + at large enough scales Q 2 . 

3. Renormalization 

We work in Feynman gauge and use dimensional regularization in D = 4 + e dimensions, applying 
the MS-scheme, if not stated otherwise. The renormalization proceeds in four steps, which we 
will briefly sketch here and refer to [^] for more details. Mass renormalization is performed in the 



on-shell scheme [10], whereas for charge renormalization we use the MS-scheme. We work in 
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an intermediate MOM-scheme for charge renormalization by requiring that the heavy quark loop 
contributions to the gluon propagator vanish for on-shell external momentum. This is necessary for 
the renormalization of the massive OMEs to cancel infrared singularities which would otherwise 
remain. Z g in this M M-scheme can be calculated using the background field method [|TI]] . Finally, 
we transform our result back to the MS-scheme for coupling constant renormalization via, rf7h, 



MOM MS 



MS Z 



+ 



nr 



(1) 



,MS~ 



with 



■>0,Q 



4 

■jT F 



-4{ -C A + C F 



32. 



P^ = -—T F C A + \5T F C F . 



(3-1) 



(3-2) 



The remaining singularities are of the ultraviolet and collinear type. The former are renormalized 
via the operator Z-factors, whereas the latter are removed via mass factorization through the transi- 
tion functions T. After coupling- and mass renormalization, the renormalized heavy flavor OMEs 
are then obtained by 

A=Z~ l Ar- 1 , (3.3) 



where quantities with a hat are unrenormalized. Note that in the singlet case operator mixing 
occurs and hence Eq. ( fO| ) should be read as a matrix equation, contrary to the NS-case. The Z- 
and T-factors can be expressed in terms of the anomalous dimensions of the twist-2 operators to 
all orders in the strong coupling constant, cf. [^, [12|] up to 0(a^). From Eq. (3.3) one can infer 
that for operator renormalization and mass factorization at 0(o? s ), the anomalous dimensions up to 
NNLO ^ together with the 1-loop massive OMEs up to 0(e 2 ) and the 2-loop massive OMEs up 
to 0(e) are needed. The 2-loop OMEs up to 0(e°) were calculated in Refs. [|, §, ||, [l3| |. Higher 
orders in e enter since they multiply Z— and T-f actors containing poles in e. This has been worked 
out in detail in Ref. [|lj], where we presented the 0(e) terms Oq , <^g S an d "Qq- ^ e terms 
sgQ an( ^ ~®&iQ were §i ven m Refs. Jl3|]. Thus all terms needed for the renormalization at 3-loops 
in the unpolarized case are known. 

Finally we would like to point out the difference between the MOM- and MS-scheme for 
coupling constant renormalization at NLO [fl3|]. Eq. ( |2.2[ ) holds only for completely inclusive 
quantities, including radiative corrections containing heavy quark loops ||. Additionally, (2.2) 
has to be applied in such a way that renormalization of the coupling constant is carried out in the 
same scheme for all quantities contributing, i.e., the MS-scheme. If one evaluates the heavy-quark 
Wilson coefficients, diagrams of the type shown in Fig. [TJ may appear as well. It contains a virtual 
heavy quark loop correction to the gluon propagator in the initial state and contributes to the terms 
L g j and H g j, respectively, depending on whether a light or heavy quark pair is produced in the 
final state. Note that in the former case, this diagram contributes to F( 2 ,l)( x i Q 2 ) m tne inclusive 
case, but is absent in the semi-inclusive <2<2-production cross section. In Refs. [|l|], the coupling 
constant was renormalized in the MOM-scheme in 0(a 2 s ) by absorbing the contributions of the 
above diagram into the coupling constant. This can be made explicit by considering the complete 
gluonic Wilson coefficient up to 0(a 2 s ), including one heavy quark, see Eq. 



C g ,2 ("/) + L gt z (tif + 1 ) + H g 2 (rif+\)=a 



MS 



(l).MS 
l Qg 
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QiQi 



Figure 1: 0(a 2 s ) virtual heavy quark corrections to H 



(2) 



4 (2), MS 



(1),MS r (l),NS 
+ A Qg S,2 



+ 1 



i+cg(«/+i; 



(3.4) 



The above equation is given in the MS-scheme. Here, the diagram shown in Fig. [I] contributes, 
corresponding exactly to the color factor Tp. Transformation to the MOM-scheme for a s , Eq. (3.1), 
yields 



C g , 2 (n/) + L gj2 (n f + 1 ) + ff gl2 («/ + 1 ) = oj 



. (l),MOM 



+ C$(n/ + 1) 



, (2), MOM 

l 2s 



+A 



(l),MOM ^(1),NS 



^r>/ + l)+<S(n/+l 



(3.5) 



In the above equation, all contributions due to diagram [l] have canceled, i.e. the color factor Tp 



does not occur at the 2-loop level in the MOM-scheme. Splitting up Eq. ( fL5| ) into H g j and L g j, 
one observes that L g j vanishes at 0(a^). The term H g j is the one calculated in Ref. [Q], which is 
the asymptotic expression of the gluonic heavy flavor Wilson coefficient as calculated exactly in 
Refs. [fi]]. It is not clear whether the same can be achieved at the 3-loop level as well, i.e., trans- 



forming the general inclusive factorization formula (2.2) in such a way that only the contributions 
due to heavy flavors in the final state remain. Therefore one should use the asymptotic expressions 
at 3 loops only for completely inclusive analyzes. This approach has also been adopted in Ref. [^] 
for the renormalization of the massive OMEs, which was performed in the MS-scheme and not in 
the MOM-scheme, as previously in Ref. [J3j] . In the NS-case a similar argument holds, which can 
be found in Ref . [g]. 



4. Calculation and Results 

The massive OMEs at 0(a 3 s ) are given by 3-loop self-energy type diagrams, which contain a 
local operator insertion. The external massless particles are on-shell. The heavy quark mass sets 
the scale and the spin of the local operator is given by the Mellin-variable N. The steps for the 
calculation are the following: We use QGRAF jl4| ] for the generation of diagrams. Approximately 
2700 diagrams contribute to all the OMEs. For the calculation of the color factors we refer to [jl^]. 
The diagrams are then genuinely given as tensor integrals. Applying a suitable projector provides 
the results for the specific Mellin moment under consideration. The diagrams are further translated 
into a form, which is suitable for the program MATAD Jl6|], through which the expansion in e is 
performed and the corresponding massive three-loop tadpole-type diagrams are calculated. We 



have implemented all these steps into a FORM-program [17] and checked our procedures against 
various complete two-loop results and certain scalar 3-loop integrals and found full agreement. 
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Applying Eq. (3.3), one can predict the pole structure of the unrenormalized results and thus 
the logarithmic terms of the renormalized OMEs. These contributions can be expressed in terms 
of the anomalous dimensions up to 3 loops, the expansion coefficients of the QCD /3-function up 
to 2 loops and the 1- and 2-loop contributions to the massive OMEs. Thus the logarithmic terms 
are known for general values of N. This is not the case for the constant term, which contains the 
genuine 3-loop contributions aj? . These are known for the fixed values of N as calculated in this 
work [JJ. 

For the OMEs A^,A^ e and A® Q the moments N = 2to 10, for A^' PS to N = 12, and for 

(3) NS (3) PS (3) 

qq Q > qq Q ' gq Q to ^ = 14 were computed. For the flavor non-singlet terms, we calculated 
as well the odd moments N = 1 to 13, corresponding to the light flavor — combinations. 1 The 
complete calculation took about 250 days of computer time. All our results agree with the predic- 
tions obtained from renormalization, providing us with a strong check on our calculation. As an 
example, we show the constant term of a^j of the unrenormalized OME Agj for N = 10 



fl e s ( 10 ) T F C A I 685850575063965696000000 81675 B " + 9075 ^ 415451499724800 ^ l"12 «"(»»> 



/ 6830363463566924692233659 563692 „ 483988, 103652031822049723, 20114890664357, 



+T F C F C A 



1 872201479486471797889957487 1286792 643396 761897167477437907 15455008277 
1 2992802509370032 1 28000000 + 81675 4 ~ 9075 ^ ~ 33236119977984000 ' 3 + 660342375 • 



,/ 247930147349635960148869654541 11808 53136 9636017 147214304991 14699237127551 
+Tf f [ 148143724213816590336000000 3025 4+ 3025 ^ 4+ 7122025709568000 ' 3 + 15689734830000 * 3 

2 I 23231 189758106199645229 123553074914173 4206955789, 



633397356480430080000 5755172290560 " 377338500 



2 / 18319931182630444611912149 502987059528463 24683221051 \ 896 , 2 / 

+ F f [~ 1410892611560158003200000 ~ 1 13048027136000 ™ + 46695639375 ~ 1485 " 3+ "/ F A l 

1505896 . 189965849 A , / 1178560772273339822317 62292104 , 49652772817, 



297277185134077151 
15532837481700000 



^ 3 + 188669950 ^ 2 I +n f^F^F\ 107649563748181000000 + 13476375 ^ 3 93391978750^- 



245025 " 188669250'-/ \ 107642563748181000000 13476375 " 93391278750 

Here Q denotes the Riemann ^-function at integer argument i and the term B4 is given by 

B 4 = -4C2ln 2 2 + ^ln 4 2-yC4 + 16Li 4 Q . (4.1) 
It appears in all OMEs we calculated and is known to arise as a genuine mass effect. 



5. Conclusions and Outlook 

We calculated all massive 3-loop OMEs for even Mellin-moments N = 2. ..10(12, 14) using 
MATAD. This confirms for the first time, in an independent calculation, the moments of the 
fermionic parts of the corresponding 3-loop anomalous dimensions Combining our results 
with this provides fixed moments of the heavy flavor Wilson coefficients of F2 in the limit 
Q 2 ^> m 2 . First phenomenological studies of the effects of our calculation are in preparation and 
will be used extending the heavy flavor treatment from 2- to 3-loop accuracy in foregoing analyzes 
@. 

'The massive OMEs for transversity have been calculated in jli|. The corresponding anomalous dimensions found 
in are confirmed in the color factors appearing in the present calculation and extended to higher values of N. 
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